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Twirl operations, which convert impure singlet states into Werner states, play an important role in 
many schemes for entanglement purification. In this paper we describe strategies for implementing 
twirl operations, with an emphasis on methods suitable for ensemble quantum information processors 
such as nuclear magnetic resonance (NMR) quantum computers. We implement our twirl operation 
on a general two-spin mixed state using liquid state NMR techniques, demonstrating that we can 
obtain the singlet Werner state with high fidelity. 

PACS numbers: 03.67.-a 



I. INTRODUCTION 

The twirl operation was introduced by Bennett et al. 
PIQ hi the context of entanglement purification of mixed 
states. The first step in many of these protocols is the 
twirl operation which converts an arbitrary mixed state 
of a two qubit system, p, into a Werner singlet state 
Pw j an incoherent mixture of e parts of a maximally 
entangled singlet state, |* _ ) = (|01) - |10>)/\/2, and 
(1 — e) parts of the maximally mixed state, 

p w (e)= E \*-)(*-\ + (l-e)j. (1) 

The key idea behind the twirl is that the singlet state 
is invariant under any bilateral unitary transformation 
of the two qubits (that is, any operation where identical 
local unitaries are applied to the two qubits), whereas 
any other state will be affected. A randomly chosen bi- 
lateral rotation should serve to average any other state, 
converting an arbitrary mixed state to a Werner singlet. 
(Note that the twirl sequence does not create new sin- 
glet, rather it preserves the singlet, while averaging out 
all other terms.) This provides a good "in principle" def- 
inition of the twirl operation, but in practice it is impor- 
tant to consider what a random bilateral rotation really 
is, and how it can actually be implemented. 

Early work in this held concentrated on reducing the 
infinite set of randomly chosen rotations down to a small 
finite set. Originally it was suggested p] that the twirl 
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could be achieved using a set of four bilateral rotations, 
although it was later shown that a set of twelve bilat- 
eral rotations are both necessary and sufficient Q . This 
approach is sensible for conventional quantum informa- 
tion processors, but with ensemble devices, such as nu- 
clear magnetic resonance (NMR) quantum computers, it 
is useful to use a different approach. In such devices 
it is in fact easier to perform rotations by angles that 
vary continuously over the ensemble than by angles from 
some small carefully chosen set. In NMR devices, for 
example, the application of a strong magnetic field gra- 
dient Q performs a rotation around the z axis with a 
rotation angle that depends strongly on the position of 
a given molecule within the spatial ensemble; integrating 
over the position within the sample is then equivalent to 
applying a random rotation. In the language of NMR 
this is often referred to as a crush gradient pulse. For 
simplicity we will sometimes describe a crush gradient as 
a random rotation around the z-axis, as the effect on the 
averaged density matrix is the same, but when consid- 
ering experimental implementations it is necessary to be 
more precise. 

The use of field gradients forms the basis of spatial av- 
eraging methods, used in NMR quantum computing to 
prepare pseudo-pure states [5|, |6j . Field gradients have 
also been used to average out error terms |7(, to project 
qubits into the Zeeman basis 8], and to simplify density 
matrices prior to partial state tomography |9(. The al- 
ternative procedure, in which one rotation chosen from 
a small set is applied, requires several experiments to be 
performed, and in the context of NMR quantum comput- 
ing is usually known as temporal averaging. 

In subsequent sections we explore different methods 
for implementing the twirl, concentrating on those which 
are best suited to ensemble techniques. We begin by 
noting that any set of bilateral rotations which performs 
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a twirl must correspond to a set of single qubit rotations 
which, when considered from the viewpoint of a single 
qubit subsystem, averages any state of a single qubit to 
the maximally mixed state. As we shall see, however, 
averaging a single qubit is a necessary but not a sufficient 
condition for a set of rotations to act as a twirl. 



II. AVERAGING A SINGLE QUBIT 

Averaging a single qubit gives some useful insight into 
the problem of twirling a two qubit state. A single qubit 
state can be represented by a ray on the Bloch sphere 
and an arbitrary rotation will move this ray over the 
surface of the sphere. Each rotation n) is defined 
by a rotation angle £ about an axis n parameterized by 
a tilt angle 9 and an azimuthal angle (/). An arbitrary 
single qubit state is completely defined by three real pa- 
rameters, which are conventionally taken as the expecta- 
tion values of the single spin Pauli operators {a x , a y , a z } 
|lfj| . For scrambling an arbitrary state to the maximally 
mixed state, we must choose a convenient set of opera- 
tions, which scrambles each of these Pauli operators indi- 
vidually. A tempting, but incorrect, approach is to apply 
random rotations about random axes, that is to take £ 
as being uniformly distributed between and 2tt, while 
the rotation axes defined by 9 and <p are uniformly dis- 
tributed over the sphere; we call this set of rotations 1Z. 

Brute force integration shows that the continuous set 
of random rotations around random axes, 1Z, does not 
completely randomize the state of a single qubit, and so 
cannot constitute the basis for a twirl operation; instead 
it reduces the Bloch vector to one third of its original 
length. In retrospect the reason for this behaviour is 
clear: random rotations around axes perpendicular to the 
original state will completely average it, while rotations 
around axes parallel or anti-parallel to the state will leave 
it unaffected. Thus the overall effect of TZ must be to scale 
down the state, rather than to average it completely. The 
significance of the scaling by one third is discussed below. 

A better definition of a random rotation is provided 
by considering the description of a rotation by means of 
its Euler angles. While many different conventions for 
describing Euler angles exist, the essential feature is that 
any rotation can be decomposed as a sequence of three 
rotations around two axes. For example any rotation can 
be achieved by the sequence of rotations 

RM R y {6) R z (£) (2) 

where rotations are applied from left to right. A random 
distribution of Euler rotations can be achieved by taking 
£ as uniformly distributed between and 27r, with the 
rotation axes defined by 9 and 4> uniformly distributed 
over the sphere; we call this set of rotations £. It might 
seem that £ is the same as TZ, but they are in fact quite 
different: in particular £ completely averages a single 
qubit, while TZ does not. This is easily seen by noting 
that the first rotation will average o~ x and a y , while the 



second and third will average a z . As any state of a single 
qubit can be written as a linear combination of these 
basic matrices and the maximally mixed state, the first 
two operations will average any state. 

It is clear from the above that a random Euler rotation 
will average any single qubit, but this process is perhaps 
excessive, as two rotations suffice. (We will see below 
that this is not true when considering a true twirl applied 
to two qubits.) The sequence of rotations 

R z (<fi)Ry(6), (3) 

where both 9 and (f> are now taken as uniformly dis- 
tributed between and 2ir, will average any single qubit; 
as before this is most simply seen by considering the re- 
sult for a x , <7 y and a z . Indeed this process can be sim- 
plified still further: as described by Bennett et al. pj 
is is possible to average a single qubit by randomly se- 
lecting from the four operations {1, a x , a y , a z }. Since a a 
is equivalent to a 180° rotation, and (neglecting global 
phases) a x — a y o z it is clear that this operation is equiv- 
alent to applying either or 180° at random, and then 
applying or 180° at random; a similar result has been 
described by Hayden et al 0] • Thus instead of choosing 
the two angles in equation|21at random from uniform dis- 
tributions, we can choose at random from two particular 
values. 

This process of replacing a continuous rotation by a 
small number of discrete values is an example of a quite 
general procedure. For rank 1 and rank 2 tensors, con- 
tinuous rotations about a fixed axis, can in fact, always 
be replaced by discrete rotations: we can "hop" between 
the angles 2irk/p, with p > 3 and k — 1,2, ... ,p, and the 
overall effect is the same as rotating continuously about 
the same axis |l2j. It is sometimes possible to use p = 2, 
but this cannot be relied on in general. 

The rotation through 2-nkjp is the member Cp of the 
cyclic group C p |13|: in this way, the continuous rota- 
tions can be replaced by discrete group operations and 
the continuous integrals over the rotation angle £, can be 
substituted by group theoretical averages, such as 

/ U& fijptfttf, ft) d!i = g U( — , n)pU\—, n). 

(4) 

The simplest discrete group that can generally be used 
for averaging is the cyclic C3 group, which is employed 
by NMR spectroscopists in the form of magic angle hop- 
ping 0, Q] , a discrete version of its continuous coun- 
terpart, magic angle spinning [T^ . The group theoretical 
framework has also been used in the context of bang-bang 
decoupling protocols [TEl IT(| , which seek to average out 
system-environment interactions by applying fast, dis- 
crete and periodic control impulses to only the system, 
with the control propagators faithfully representing the 
members of a discrete cyclic group C n . 

This approach hints at another method for performing 
complete averaging of a single qubit. We have already 
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noted that random rotations around a random axis do 
not completely average a single qubit, but simply rescale 
its length by one third. This process is, as usual, equiva- 
lent to applying a rotation through an angle of 0°, 120° 
or 240° around a random axis. Clearly a 0° rotation can 
have no effect, while the effect of the 120° and 240° rota- 
tions must be the same. From this it can be deduced that 
a 120° rotation around a random axis will completely av- 
erage a single qubit, and this is indeed the case. 

Finally we turn to the issue of practical experimental 
implementations. For a conventional quantum informa- 
tion processor the obvious approach is to use the smallest 
discrete set of operations; for the case of averaging a sin- 
gle qubit this is the set of four rotations described by 
Bennett et al. With an ensemble processor, such as an 
NMR quantum computer, it is better to use a procedure 
based on continuous rotations, corresponding to spatial 
averaging. Clearly the best approach is to apply random 
rotations around two orthogonal axes, such as z and y. 
A strong magnetic field gradient, denoted G, will effect a 
random rotation around the z axis, and rotations around 
other axes can be achieved by combining gradients with 
single qubit gates (radiofrequency pulses) which can be 
treated as rotating the axis system. Thus the sequence 

G90.G (5) 

will completely average a single spin. Note that this pro- 
cedure is not completely equivalent to equation |3 as it 
should strictly speaking be followed by a 90_ x pulse to 
rotate the axes back to their original positions, but as 
the maximally mixed state will not be affected by the 
rotation this final stage may be safely omitted. 

The discussion above assumes that the gradient pulses 
are instantaneous, so that their sole effect is average the 
spin state and evolution under any background Hamilto- 
nian can be ignored. In fact, gradient pulses take a finite 
time, but as the gradient and background Hamiltonians 
commute they can be treated as an instantaneous pulse 
followed by a finite period of evolution under the back- 
ground Hamiltonian. For a single NMR spin this is not 
a problem, as the maximally mixed state does not evolve 
under the background Hamiltonian; for a two qubit sys- 
tem, however, it is necessary to consider this problem 
more carefully. 

III. TWIRLING TWO QUBITS 

We now turn to the problem of implementing a full 
twirl on two qubits. This is clearly related to the prob- 
lem of averaging a single qubit, but is more complicated. 
In particular, any twirl procedure will average a single 
qubit, but not every averaging procedure can be con- 
verted to a twirl. For example, the sequence of four op- 
erations {l,a x ,a y ,a z } suggested by Bennett et al. [I| 
does average a single qubit, but does not effect a twirl; 
instead it is necessary to use a set of 12 operations Q 
corresponding to the rotational symmetry elements of the 



tetrahedral group T [l^. The group theoretic al i ustifi- 
cation for using these rotations is discussed in |17| . 

Turning to the other averaging methods discussed 
above, only the method based on random Euler rotations 
(£ ) will correctly implement a twirl; the other methods 
will only partially average the state. In particular the se- 
quence of two random rotations around orthogonal axes 
has already been studied |9j, and shown to implement a 
partial twirl, in which the density matrix is reduced to 
Bell diagonal form, but not to a Werner state. 

For an ensemble quantum information processor, it 
might seem that the best strategy for performing a full 
twirl would be to perform random Euler rotations. Un- 
fortunately this approach cannot be implemented using 
gradients or similar methods, as these generate rotations 
with a uniform distribution (between and 27r) for each 
rotation angle, while the Euler rotation method requires 
that the angles 9 and <p be uniformly distributed over a 
sphere. Thus it might seem that twirl operations can- 
not be easily implemented by ensemble methods. It is, 
however, possible to do this by using sequential random 
rotations around three different axes. 

This procedure is a development of the partial twirl 
described above. Consider again the sequence of rota- 
tions in equation [5j These two rotations constitute a 
partial twirl 0, and leave the density matrix Bell di- 
agonal with equal populations of |<i> + ) and |3>~), where 
|$ ± ) = (| 00) ± |ll))/\/2. The role of the third rotation 
is to scramble these states with the sole remaining un- 
desired term, |* + ) = (|01) + |10))/\/2, leaving a Werner 
singlet state. This can be achieved by rotating about an 
axis at an angle tp — arccos(l/v / 3) ~ 54.74°, comm only 
known in NMR studies as the "magic ang le" \muM- 
Our final twirl sequence is then 

90^6 54.74^0. (6) 

As was the case for the partial twirl [l| it is necessary 
to choose the lengths of the gradient pulses as multiples 
of 1/6, where 6 is the difference between the Larmor fre- 
quencies of the two spins, to refocus evolution under the 
background Hamiltonian, which occurs at this difference 
frequency. 

If desired this twirl can be implemented using discrete 
steps by replacing each random rotation with a three step 
averaging procedure, 

{Z?}90 x {Z2}54.74 x {Zl}, (7) 

where {Z?} indicates a set of z rotations with rotation 
angles 2ir a/b, where a = 0, 1, ... b— 1. This discrete twirl 
requires twenty-seven steps, but the number of steps can 
be reduced to eighteen, using the sequences 

{Z?}90 X {Z$} 54.74° {Z%}, (8) 

or 

{Z?}W x {Z%}M.n x {Zl}. (9) 
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It is not possible to reduce the number of steps further 
while using sequential rotations around three axes. 

It is instructive to compare this discrete twirl with the 
twelve step version of Bennett et al. @ . Their twelve ro- 
tations include the identity operation, ir rotations about 
the principal x, y and z axes and rotations by ±27r/3 
around the four body diagonals, each of which is at the 
magic angle with the cardinal axes. (These are the twelve 
rotations that leave a tetrahedron invariant.) Clearly our 
rotations are closely related to theirs; the use of sequen- 
tial rotations means that the number of steps required 
is larger, but the sequential approach may be simpler to 
implement in practice. 

IV. AN NMR EXPERIMENT 

Finally, we demonstrate the experimental implemen- 
tation of our ensemble twirl operation as described in 
equation El Our implementation is based on a homonu- 
clear system, comprising the two 1 H nuclei of cytosine 
dissolved in D2O; we use the product operator descrip- 
tion [13 for the NMR states and pulse sequences, and 
label the two spins as I and S. The Hamiltonian of our 
system is then 

H/fr=2iTViI z + 2irv s S z +irJ2I z S z (10) 

where i>j and vg are the resonance offset frequencies of 
spins / and S, and J is the IS spin-spin coupling con- 
stant (assuming weak coupling), all measured in Hertz. 
NMR experiments were performed on a Varian Unity 
INOVA 600 MHz spectrometer with the spectrometer 
frequency placed between the two cytosine resonances; 
for our system vj = 457.9 Hz, v$ — —457.9 Hz and 
J = 7.2 Hz. 

The experiment began with the preparation of a 
generic mixed state of the two qubits containing a con- 
tribution from the singlet state as well as a wide range of 
other terms. The preparation sequence (which was ap- 
plied to the thermal equilibrium state, I z + S z ) involves 
shaped pulses selectively exciting the spins and separated 
by a delay r, 

A = 60/j, - t - 30S y . (11) 

Shaped pulses were implemented using strongly modu- 
lated composite pulses as described by Fortunato et al. 
|19| . The deviation density matrix of the prepared state 
p has a singlet component proportional to 

— sin (ttJt) sin (2ttl>jt) . (12) 
8 

Clearly the amount of singlet in p can be easily controlled 
through choice of the delay r. Note that as we are using 
a deviation density matrix description of a highly mixed 
state it is possible for the singlet component to be nega- 
tive. 



TABLE I: Pulse sequences for the first set of experiments, 
which demonstrate the effects of the three stages of our twirl 
operation. See main text for details. 



Stage 


Pulse sequence 









A - Acq 
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A — Gi — Acq 






2 


A - Gi - 90* - 


G 2 


- Acq 


3 


A - Gi - 90* - 


G 2 


- 54.74° - G 3 - Acq 



We performed two sets of experiments on states drawn 
from this family of states. The first set uses a fixed initial 
state and demonstrates a stepwise progression through 
the twirl sequence, with the resulting state becoming 
closer to a Werner state at each stage. The second set 
uses a range of initial states, with different singlet frac- 
tions, and shows that the twirl sequence works equally 
well over this range. 

The first set of experiments, listed in table HJ com- 
prises a stepwise progression through the twirl sequence, 
equation Each experiment uses the same initial state, 
with r = 69.3 ms, a value chosen to maximize the singlet 
fraction. The three stages of the twirl correspond to the 
application of three different crush gradients, and can be 
characterized as follows: stage 0, do nothing; stage 1, 
apply a single crush gradient; stage 2, apply the partial 
twirl sequence [j|; stage 3, apply the full twirl sequence. 
Two different measurements are performed to character- 
ize the state after each step: simple direct acquisition 
(which should show no signals for a Werner singlet state) 
and acquisition after a selective excitation pulse, which 
converts the singlet state into NMR observable terms. 
The selective pulse was implemented using the sequence 

90 45 - [^] - 90 180 , (13) 

where 5 = 915.8 Hz is the frequency separation between 
the two resonances as before. This sequence is based on 
iump-and-return sequences described in previous work 
[1,|2(J], an d when applied to a singlet state, results in the 
observable NMR terms, 

~(-2I x S z + 2I z S x ), (14) 

corresponding to a pair of antiphase doublets with equal 
and opposite intensities. The resulting spectra are shown 
in Figure ^ The left hand column of this figure shows 
spectra from direct acquisition and the right hand column 
shows spectra obtained with the excitation pulse followed 
by acquisition. The four rows correspond to observation 
after 0, 1, 2 or all 3 stages of the twirl sequence. 

The initial state p contains many different components, 
and the observed spectra (top row) are complicated, re- 
flecting this fact. After the first stage of the twirl se- 
quence (the first field gradient) all components which 
are directly observable by NMR are crushed (averaged to 
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Stage 



Direct detection 



Detection after Excitation 



1 



r 



i 



FIG. 1: Experimental spectra depicting an implementation of an ensemble twirl. Spectra were obtained using the pulse 
sequences listed in table|l] The four rows correspond to spectra acquired after 0, 1, 2 or 3 stages of the twirl sequence, while the 
left and right columns correspond to direct signal acquisition and acquisition after the selective excitation sequence. The ideal 
result at the end of the full twirl (last row) is no signal with direct detection (left column), and a pair of antiphase doublets 
with equal and opposite intensities after a selective excitation sequence (right column). 



zero), and so no signal is visible in the direct detection 
spectrum, but many other components remain, indicated 
by the variety of signals seen after excitation. 

The second stage of the twirl (the 9Q X pulse and the 
second gradient) removes most of these terms, produc- 
ing a Bell diagonal state with equal populations of |$ + ) 
and \&~). As before there is no signal in the left hand 
spectrum, while the spectrum on the right contains two 
antiphase doublets with clearly different intensities. This 
intensity difference arises from the imbalance between the 
|*+) state and the two 1$*) states Hi- 
After the third and final stage of the twirl (the 54. 74^ 
pulse and the last gradient) these two antiphase doublets 
have equal and opposite intensities, characteristic of a 
Werner singlet state. The final intensity of these dou- 
blets, compared with a calibration spectrum (not shown) 
acquired from the thermal state, is consistent with the 
twirl preserving the fraction of the singlet state as ex- 
pected. 

This last point is explored in more detail in the second 
experiment, in which the full twirl is applied to a range of 
initial states with different amounts of the singlet state. 
The size of the singlet component for our family of states 
is given in Eauation ll2l and shows two kinds of sinusoidal 
modulation with the variable delay r: a fast variation, 
arising from the offset frequency vj, on top of a slow 
variation due to the coupling J between the two qubits. 
If we choose r to be close to 1/(2 J), we are near the 
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FIG. 2: Spectra from Werner singlet states obtained by 
twirling impure states with varying singlet fractions; only 
one component of one antiphase doublet is shown. Successive 
spectra are obtained by incrementing the variable delay in the 
preparation sequence, Eg. 1111 in steps of 1/(10^/) = 218 (is, 
which causes the fraction of singlet in the initial state to be 
modulated, with a period of ten steps. 



maximum of the slow J modulation, and the effect of 
varying r is dominated by a sinusoidal variation arising 
from uj. 

In this way we can produce a range of density matri- 
ces p with varying amounts of singlet, together with other 
terms, each of which can be twirled to produce a Werner 
state. This state can then be observed using a selective 
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excitation pulse prior to acquisition. In each case the ex- 
pected spectrum is a pair of antiphase doublets, with the 
intensity of the signal showing a sinusoidal modulation 
at the frequency z/j. The results of this experiment are 
shown in Fig. [2J which depicts the intensity variation in 
the righthand component of the lefthand doublet as r is 
varied around a value of 1/(2 J) = 69.3 ms, with an in- 
crement between successive spectra of l/(10Vj) = 218 /is; 
equivalent effects can be seen for the other three compo- 
nents of the NMR signal. As expected a sinusoidal mod- 
ulation of the signal is clearly seen, and the observed 
modulation period of ten spectra is exactly as expected. 
The small out of phase signals observable near the zero- 
crossings of the sine wave can be ascribed to the effects 
of spin-lattice relaxation during the gradient pulses. 



turn computer. With a conventional quantum computer 
the implementation requiring the smallest number of dif- 
ferent bilateral operations is the set of 12 rotations pre- 
viously described 0, but our 18 step and 27 step aver- 
aging procedures require a smaller number of elementary 
operations and may be simpler to implement in practice. 
With an ensemble quantum computer, such as an NMR 
device, it can be simpler to replace the discrete averaging 
procedure by continuous averaging, exploiting the ensem- 
ble nature of the system. We have developed a scheme, 
involving the application of three successive crush gra- 
dients separated by RF pulses, which is well suited to 
NMR quantum computers, and have demonstrated that 
its experimental performance is consistent with theoret- 
ical expectations. 
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